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discussed.  A  model  of  mixing  by  wind-generated  wave  motion  is  developed 
and  is  tested  along  with  two  other  models;  and  the  results  are  discussed. 
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I.  Introduction 

The  plienoinenon  jf  mixing  of  the  surface  c/f  r^uch  of  the 

world’s  ocean  area  has  been  studied  for  many  years.  Increasing 
naval  o^^erational  interest  in  the  depth  of  the  mixed  layer  nov;  demands 
that  the  groundwork  established  by  many  oceanographic  investigators  be 
built  upon,  and  that  intensive  effort  be  directed  toward  the  goal  -  a 
usable  forecast  of  mixed-layer  depth  over  broad  expanses  of  the  seas. 

Mixed-layer  depth  is  the  depth  below  the  water  surface  to  xjhich 
mixing  fi^s  established  an  essentially  isothermal  temperature  distri¬ 
bution.  The  lower  boundary  of  the  mixed  layer  is  the  thermocline, 
normally  a  thin  layer  or  interface  of  large,  negative  vertical  temper¬ 
ature  gradient. 

Large  variations  of  mixed-layer  depth  occur  in  space  and  time. 
Observation  indicates  that  a  full  spectrum  of  time  variation  ranging 
from  an  annual  cycle  to  short  term  fluctuations  of  a  few  minutes 
duration  is  to  be  expected  at  a  fixed  geographic  location.  The  annual 
cycle  is  at  present  the  most  regular  fluctuation  of  the  mixed-layer 
depth  to  be  identified^  and  is  in  general  the  largest  in  magnitude. 

The  annual  trend  of  the  mixed  layer  is  a  sinusoidal  wave  which  appears 
to  be  most  closely  related  to  the  heat  balance  of  the  surface  water 
layers,  i.e.  mixed-layer  depth  is  inversely  proportional  to  the  sta¬ 
bility  of  tne  surface  water  layers.  This  is  not  to  say  that  other 
periodic  fluctuations  of  mixed- layer  depth  such  as  lunar- tidal  or 
diurnal  heating  and  cooling  are  not  present,  but  paucity  of  data  and 
the  nxignitude  of  apparently  random  fluctuations  makes  their  identifica¬ 
tion  difficult. 

instantaneous  values  ol  mixed-layer  depth  may  be  looked  upon  as 
the  result  of  various  perturbations  superimposed  upon  the  mean  annual 
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wave.  The  most  si^iiifixant  perturbations  appear  to  be  the  integrated 
cfrcct  of  meteorological  factors  operating  to  alter  tliu  stability  of 
Lac.  surface  v^atcr  layers.  \/iiids,  evaporative  cooling,  albedo,  moisture 
coiitent  of  the  air,  and  insolation  are  a  lev?  of  these  meteorological 
factors.  Added  to  the  effects  of  the  meterological  perturbations  are 
the  a^.pareiitly  random  motions  of  internal  vjaves  along  the  thermocline. 

Geographical  variations  of  mixed*- Layer  depth  are  associated  vjith 
the  ba th;;’mctry  of  the  ocean  basins,  the  advective  and  mixing  effects 
of  tlie  permanent  currents  and  ocean  tides,  and  the  climatological 
regimx^s.  The  enorraous  problems  of  data  collection  and  observational 
density  militate  against  anything  but  qualitative  estimates  of  the 
effects  of  space  variation  upon  luired-layer  depth  at  the  present  time. 

The  purpose  of  this  paper  is  to  evaluate  the  raetcorological  factor 
of  wind  in  its  contribution  to  mixing  of  tlie  surface  waters  at  ocean 
station  Papa  (50li  -  145\0  for  the  year  1953. 
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2. 


Tiic  meciic.iisin  of  ’./mcl- in<'luccd 


Cell  cCL  e  rectly  to  ,  roducc  r.ioineutuin  troasfei:  across  tu*^ 
air-sea  i.it«_cfacc  by  means  of  wind-driven  wave  liiotiou  or  by  drift 
currents  ClZI  •  fs  postulated  that  a  given  e'lmount  oi  v;ind  stress 

will  produce  raouioii  of  the  surface  layers  of  \?ater,  giving  rise  to 
lair.lug  of  the  v;ater  layers  until  stability  j.ii  the  rcgi.jn  of  the  ther- 
mocllne  dampens  particle  notion  to  such  an  extent  that  equilibrium 
between  th^  mi:c-.ng  forces  and  tne  stabilizing  forces  is  reached,  and 
/crtical  lulriug  can  proceed  no  further.  The  aLX)unt  of  stress  exerted 
by  tae  \;iad  upon  tlie  sea  surface  is  kno’./n  to  depend  u^>o'n  tne  speed  of 
the  wind,  the  stability  of  the  air,  and  the  nature  of  the  air-sea 
interface  [^2^,  The  stability  of  the  ther.Qocline  (the  limit  of  mixing) 
depends  upon  the  V7iiif]  stress  C3]3  >  ihe  heat  balance  in  the  mixed 

layer . 

Convective  mixing  is  the  process  by  v/hich  evaporation  and  loss 
of  heat  to  the  atmosphere  causes  the  surface  v/ater  layers  to  become 
unstable,  resulting  in  overturning  and  mixing.  Cind  speed  is  a  factor 
in  the  flux  of  water  vapor  and  sensible  heat  across  the  a>.r-sea  inter¬ 
face.  Convective  nixing  is,  in  general,  most  efficient  when  heat 
losses  from  the  surface  mixed  layer  are  greatest. 

In  order  to  isolate,  Insofar  as  possible,  the  effect  of  V7ind- 
driven  mechanical  mixing  for  tne  purpose  of  study  in  tals  section, 
data  from  the  months  of  maximum  sea- surface  neating  were  chose;,  from 
the  records  of  the  pacific  Oceanic  Group  for  ocean  station  Pa^^a,  1553, 
as  ex.aiu^lcs  of  tnis  process. 


3.  x.ic  r^la tioncliii^  of  wlad  S;>aod  Lo  nixcd-laycr  n^pta 

i\ic  data  \}ci'^  ov*j,ariiiiCcl  lato  time  se^'ics  of  laixcc -laye^*  c\.[  Lit.  end 
\viac  s^ced,  and  into  scatter  diaj;rar.is  C'f  tiic,  sere  variables.  Liacc 
\.ave  motion  depends  on  duration  of  wind  rather  than  instantaneous  values, 
the  mean  of  the  eight  available  daily  wind  observations  V7as  taken  as  the 
representative  wind  for  each  day.  The  selection  of  a  representative 
mixed- layer  depth  tor  a  day  was  a  greater  problem.  In  general,  only 
tv7o  bathythermograph  observations,  one  at  0200  and  the  other  at  17u0 
Grceiwich  time,  were  available.  In  addition,  as  v^7ill  be  discussed 
later,  considerable  fluctuation  of  MLD  about  the  daily  mean  value  is 
observed,  and  there  is  no  certainty  that  the  mean  of  the  0200  and  1700 
observations  would  be  truly  representative  of  the  mean  for  the  day. 
Therefore,  the  maximum  observed  MLD  for  the  day  v/as  used.  MLD  v;as 
taken  to  be  the  depth  below  the  sea  surface  at  whech  the  temperature 
trace  of  the  water  ceased  to  be  isothermal;  in  general  this  is  the  top 
of  the  thermocline .  This  definition  permits  a  zero  value  for  MLD. 

Figure  1  is  a  time  series  of  mean  daily  wind  speed  and  mriximum 
observed  daily  MLD  at  ocean  station  Papa  for  the  month  of  June,  1958. 

The  two  curves  are  very  similar,  increasing  and  decreasing  together. 

Close  inspection  reveals  that  the  MLD  curve  lags  behind  the  wind  speed 
curve  by  about  24  hours,  and  that  most  fluctuations  of  wind  speed  are 
follov7ed  the  next  day  by  a  similar  fluctuation  of  MLD.  Another  interest¬ 
ing  characteristic  of  the  two  curves  is  that  not  only  does  MLD  increase 
V7ith  increasing  v/ind  speed,  it  also  tends  toward  zero  as  wind  speed 
decreases . 

Figure  2  is  a  scatter  diagram  obtained  by  plotting  maximum  daily 
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MLD  against  the  previous  day*s  mean  wind  speed  for  the  months  of  June, 
July,  August,  and  September,  1958.  Correlation  coefficients  for  the 
four  months  were  found  to  be  0.94,  0.76,  0.41,  and  0.68,  respectively. 
Disregarding  for  the  moment  the  August  points,  it  is  seen  that  the 
remaining  points  fall  into  two  groups;  those  of  June  and  July;  and 
those  of  September.  The  August  points  can  be  divided  into  two  groups 
also;  those  associated  with  the  June  and  July  group,  and  those  assoc¬ 
iated  with  the  September  group.  The  separation  of  the  points  into  two 
groups  is  in  accord  with  the  two  different  thermal  regimes  which  are 
represented.  The  June  and  July  points  represent  a  period  where  the 
water  is  being  actively  heated,  while  the  September  points  represent 
a  period  of  beginning  heat  loss  by  the  water;  the  Augi4st  points  are 
transitional  between  the  two  periods.  The  basis  for  this  division  is 
that  the  observed  sea  surface  temperature  reaches  its  maximum  for  the 
year  in  late  July  and  begins  its  annual  cooling  trend  in  mid-August. 
This  difference  in  regimes  may  explain  in  part  why  the  correlation 
coefficients  for  June  and  July  are  so  much  larger  than  the  August  and 
September  coefficients.  In  all  four  of  the  months  increasing  winds 
drive  the  mixed  layer  deeper,  but  only  under  conditions  where  excess 
heat  is  being  supplied  is  there  a  tendency  for  the  mixed  layer  to 
approach  zero  under  decreasing  wind  conditions. 

Linear  regression  equations  were  calculated  for  these  scatter 
diagrams  and  found  to  be: 

MLD  =  1.14  V  (June)  (MLD  iu.  meters,  V  in  knots) 

MLD  =  0.94  V  1.2  (July) 

MLD  =  1.07  V  26.1,  (Sept) 
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respectively.  The  small  difference  between  the  June  and  July  slopes 
is  attributed  to  the  slightly  more  stable  condition  of  the  July  water 
when  heating  was  at  a  maximum. 

The  fact  that  mixing  of  the  surface  water  can  depend  not  only 
upon  wind  speed,  but  upon  convection  as  well,  is  illustrated  in  figure 
3.  This  figure  is  a  scatter  diagram  of  the  same  variables  shown  in 
figure  2,  but  for  the  month  of  November,  1958,  Here  a  correlation 
coefficient  of  -0.58  is  obtained.  In  November  there  is  rapid  cooling 
of  the  surface  water,  the  temperature  decreasing  some  four  and  one- 
half  degrees  fahrenheit  in  two  weeks.  In  this  autumn  month,  clearly, 
convection  has  completely  overridden  the  effects  of  wind  mixing  alone, 
and  is  independent  of  wind  speed. 
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A.  Time  lag  of  the  mixed  layer  produced  by  a  given  wind 

As  pointed  out  in  section  3,  changes  of  MLD  lag  changes  of  wind 
speed.  Figure  A  is  a  time  series  of  wind  and  MLD  observations  taken 
at  ocean  station  Papa  on  16  to  19  June,  1958.  MLD  observations 
were  taken  hourly,  while  wind  observations  were  taken  at  three-hour 
intervals.  Both  curves  were  smoothed  by  three-hour  overlapping  sums, 
and  the  MLD  curve  is  plotted  with  an  18-hour  lag  from  the  wind  speed 
observations  to  obtain  good  agreement  of  the  fluctuations  on  both 
curves.  Actually  the  best  agreement  would  have  been  obtained  by  a 
12-hour  lag  for  the  shallow  MLD  (smaller  wind  speed)  with  a  larger  lag 
for  the  deeper  MLD  values.  A  lag  between  observed  \7ind  speed  and  the 
resulting  mixing  of  the  surface  water  is  in  accord  with  the  theory  of 
wave  generation,  v/hich  requires  a  minimum  duration  time  for  the  wind 
to  have  blown  in  order  to  produce  waves  of  a  certain  size.  In  general, 
the  required  duration  time  of  the  v/ind  increases  with  increasing  size 
of  V7ave  produced  [j^  ;  by  analogy,  the  lag  of  MLD  from  wind  should 
increase  with  increasing  wind  and  mixed-layer  depth. 
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5. 


ivandom  fluctuations  of  MLD 


The  depth  of  the  mixed  layer  was  observed  to  var'y  several  meters 
in  a  random  fashion  between  hourly  measurements.  The  laagnitude  of  the 
fluctuations  could  not  be  made  to  correlate  with  any  observed  meteoro¬ 
logical  parameter,  but  the  fluctuations  did  appear  to  increase  in  size 
\7ith  depth. 

Figure  5  shows  the  mean  hourly  fluctuation  of  ten  series  of 
observations  plotted  as  a  function  of  the  mean  mixed-layer  depth  of 
each  scries.  The  individual  series  contained  from  24  to  48  bathythermo¬ 
graph  observations,  and  were  spaced  at  approximately  monthly  intervals 
over  the  year. 

Schule  [3Z] others  have  observed  similar  short-term  fluctuations 
of  MLD  at  other  localities  and  conclude  tnat  they  are  caused  by  internal 
v;avcs  along  the  thermocline. 

The  importance  of  the  random  fluctuations  is  tv/o-fold.  A  single 
measurement  of  MLU  will  not  necessarily  be  representative  of  the  mean 
liLD  for  any  given  day.  A  series  of  observations  should  be  taken  each 
day  v;hcn  verifying  computed  values  of  raixed-laycr  depth.  The  other 
important  aspect  of  internal  wave  motion  along  the  thermocline  is  that 
internal  waves  may  provide  one  mechanism  of  vertical  mixing.  Ball  £^6 3) 
describes  an  experiment  in  which  internal  waves  were  induced  along  the 
density  discontinuity  betv;een  two  water  layers.  The  waves  were  observed 
to  become  sharp-crested  and  filaments  of  the  denser  water  were  dravm 
out  from  the  sharp  \7ave  crests  and  mixed  into  the  upper  layer.  The 
reverse,  mixing  dovaiwards,  was  not  observed  tc  occur. 
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6. 


Theoretical  models  of  wind  mixing 


At  the  start  of  the  present  investigation  a  search  of  the  literature 
produced  only  one  theoretical  model  to  explain  the  mechanism  of  wind 
mixing.  Munk  and  Anderson  OU  developed  this  model  incorporating  the 
concept  of  vertical  Austausch  coefficients  in  the  ocean,  solving  a  system 
of  five  equations  to  obtain  the  relationship  of  wind  speed  to  the  depth 
of  the  mixed  layer;  they  obtained  results  to  within  an  order  of  magni¬ 
tude  of  observed  data,  and  concluded  that  convective  mixing  was  at 
least  as  important  as  v;ind  mixing.  This  is  not  surprising  in  view  of 
Laevastu*s  observation  Ql']J  that  the  vertical  Austausch  coefficient 
varies  by  more  than  three  orders  of  magnitude  with  space  and  time  in 
the  oceans  and  is  at  present  neither  measureable  nor  predictable. 

The  author  approached  the  problem  of  wind  mixing  from  the  stand¬ 
point  of  particle  motion  produced  by  wind-driven  v/aves.  Essentially, 

V7hat  is  required  is  a  way  to  link  wind  speed  to  the  motion  of  v/ater 
particles  in  the  surface  water  layers  resulting  from  wind-produced 
v;aves.  having  accomplished  this,  a  limit  remains  to  be  set  on  the 
effectiveness  of  the  particle  motion  to  produce  mixing.  The  relation¬ 
ship  of  \7ind  speed  to  the  various  surface  wave  parameters  (significant 
wave  height,  period,  and  length)  was  obtained  using  Neumann’s  spectrum 
J ,  snd  the  subsurface  particle  motion  was  in  turn  calculated  from 
the  theory  of  simple  Airy  waves.  The  maximum  limit  of  mixing  due  to 
particle  motion  was  assumed  to  be  the  depth  were  the  buoyant  force  of 
the  denser  water  below  the  thermocline  exactly  opposed  the  dov7nv7ard 
force  of  the  denser  particles  in  circular  orbit.  An  equation  (hereafter 
referred  to  as  Geary’s  model)  of  the  follov7ing  form  was  derived: 


The  depth  at  which  equilibrium  will  be  established  depends  upon  the 
stability  of  the  tnerniocline  (the  density  difference  across  the 
thermocline)  and  upon  the  characteristics  of  the  surface  wave  produci  g 
the  motion.  The  necessary  wave  characteristics  can  in  turn  be  obtained 
as  a  function  of  wind  speed  alone,  assuming  fully-developed  seas. 

During  the  present  investigation  the  author  became  aware  of  the  fact 
that  Laevastu  had  approached  the  problem  of  wind  mixing  from  a 

similar  viewpoint.  The  essential  difference  betv/een  the  two  models  is 
that  Laevastu  set  an  arbitrary  limit  to  the  effectiveness  of  subsurface 
particle  motion  in  accomplishing  mixing.  dis  limit  x/as  the  depth  ^niere 
the  diameter  of  the  particle  orbits  was  tea  centimeters  or  less. 

Neumann  (373  >  an  extension  of  his  surface  wave  spectrum  to 
subsurface  motion,  derived  a  relationship  giving  the  ratio  of  potential 
energy  at  some  depth  to  that  of  the  surface  waves,  and  proposed  that  at 
the  depth  where  the  ratio  was  five  percent  or  less,  subsurface  particle 
motion  could  be  considered  negligible.  Neumann  did  not  apply  this 
concept  to  the  problem  of  mixed -layer  depth,  but  the  l^vel  of  v;ave 
energy  present  at  any  depth  must  have  some  influence  on  the  v7ork  done  in 
mixing  at  that  depth.  Consequently,  Neumann* s  relationship  was  evaluated 
along  \;ith  the  models  of  Laevastu  and  Gear}'. 
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7. 


Fit  of  models  to  observed  data 

Tiie  appendices  coutaiu  evaluatioas  of  the  thieo  mixing  eiodcls  for 
various  \;iud  speeds.  Curves  1,  2,  3,  and  A  of  figure  6  are  plots  of 
tables  1,  2,  3,  and  4  respectively,  from  tae  appendices.  The  stippled 
area  represents  the  field  ol  scatter  ot  the  observed  data  for  the 
months  oi  June,  July  and  the  first  halt  of  xiugust,  1958,  some  76  obser¬ 
vations  in  all. 

Curve  2  is  tne  Lcevastu  equation  for  HLr  using  his  equation 
0->  ??•  70^  for  computing  significant  \/ave  height.  it  can  be  seen 

chat  curve  2  forms  an  upper  limit  on  tne  observed  scatter  of  points. 

Curve  3  is  the  Laevastu  equation  for  MLl),  but  using  (V]  for  the  calcu¬ 
lation  of  significant  wave  height.  It  is  evident  that  the  method  of 
calculating  surfac«^  wave  characteristics  has  an  importc' at  effect  upon 
tne  values  obtained  for  mixed-layer  depth. 

Curve  1  is  Geary’s  equation,  and  it  Is  in  surprisingly  good 
agreement  with  curve  3.  An  explanation  for  this  agreement  cannot  be 
suggested  at  the  present  time. 

Curve  4  is  the  MLD  equation  based  on  Neumann’s  ratio  of  potential 
energies,  ^:nd  it  can  be  seen  that  this  curve  fits  the  observed  data 
better  than  the  others.  It  may  well  be  that  a  calculation  of  energ}^  at 
depth  is  a  better  measure  of  mixing  efficiency  than  are  approaches  using 
particle  acceleration  or  the  geometry  of  motion. 

Of  some  interest  is  the  fact  that  all  of  the  curves  deviate  more 
from  the  trend  of  the  observed  scatter  points  at  lowc::  wind  speeds, 
which  may  indicate  deficiencies  at  low  wind  speeds  in  the  wave  generation 
equations.  Munk  and  others  have  proposed  a  ’’critical  vjind  speed”  (about 
13  knots),  above  and  below  which  wave  generation  by  wind  proceeds 
differently. 
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0.  Variation  of  observed  data  frou  computed  values 

Factors  operating  to  cause  variation  of  observed  i;iad  and  MLD 
data  frora  theoretical  and  empirical  relationships  arc  of  three  kinds. 

First  are  tUose  that  have  to  do  with  data  collection  and  processing. 
ijath>  thermograph  observations  at  ocean  station  Papa  were  not  actually 
taken  at  a  fixed  location,  but  over  an  area  of  many  square  miles.  The 
basic  grid  is  ten  nautical  miles  on  a  side,  and  observations  are 
reported  fiwm  several  different  grid  positioiiS  in  the  course  of  a  month. 
The  best  accuracy  to  which  the  bathythermograph  lata  could  be  read  is 
on  the  order  of  ~  o.^e  meter.  It  is  entirely  possible  that  the  differ¬ 
ences  in  geographical  location  at  v/nich  the  observations  were  taken 
could  account  tor  the  scatter  of  observed  data  from  eomputed  values. 

A  second  factor,  concerning  only  the  mixing  models,  stems  from 
over s5nuplif ication  of  the  vertical  i.iixiug  process  and  from  imperfections 
in  wind-produced  \7ave  theory.  The  author’s  assumption  (appendix  I) 
of  a  tco-layer  system  is  an  oversimplification  of  actually  existing 
density  gradients.  That  systematic  part  of  the  variation  of  the  mix¬ 
ing  models  from  observed  data  may  be  ascribed  to:  1)  wind  may  not  build 
\7aves  exactly  according  to  the  theory  used:  2)  use  of  a  24-hour  mean 
vrind  instead  of  an  integrated  vjind  could  produce  a  systematic  error  of 
the  type  observed  betv/een  the  scatter  points  and  curves  on  figure  4. 

The  third  factor  is  internal  wave  motion.  Random  hourly  variation 
of  MLD  \7as  found  to  exist  (section  5)  and  was  calculated  to  be  an  increas¬ 
ing  function  of  mixed-layer  depth.  Mean  hourly  variation  of  liLD  increased 
from  nearly  three  meters  at  shallow  depths  to  «aboat  five  or  si::  meters 
at  large  MLD  values.  Also,  the  standard  deviation  of  MLD  from  the 


monthly  mean  \7as  found  to  be  a  function  of  \.ind  speed,  depth,  and 
hert  balance.  It  is  seen  that  internal  wve  motion,  in  addition  to 
observational  factors,  could  account  for  all  of  the  random  scatter  of 
the  plotted  points. 


13 


9. 


Conclusions  and  recomiaendations  lor  future  research 


Wind  is  an  important  factor  In  determining  the  depth  to  v:hich  the 
surface  v/aters  arc  mi::cd.  The  mechaaisra  of  udi^iing  by  wincl-drivcu  wave 
motion  must  e::ist  throughout  the  seasons,  but  is  dominant  only  \;hen  the 
miacd-layer  depth  tends  toward  zero  in  the  absence  of  mixing.  In  gen¬ 
eral,  this  coi.dition  will  occur  during  the  surijmer  montlis  \;hcn  heating 
of  the  surface  water  is  maximum,  and  convective  mixing  is  at  a  minimum. 

Equations  incorporating  the  subsurface  particle  motior.  of  vdnd 
\;aves  cm  be  called  upon  to  hindcast  much  of  the  observed  daily  and 
longer-iaterval  changes  in  mixed-layer  depth  at  ocean  station  Papa 
during  the  smiuncr  of  1958.  Although  turbulent  mixing  due  to  the  vertical 
velocity  shear  of  drift  currents  also  must  take  place,  this  effect  seems 
to  be  masked  by  wave  mixing,  since  the  equations  for  wave  raixrng 
generally  give  values  of  MLD  greater  than  those  actually  observed  at 
wind  speeds  greater  than  about  12  knots.  it  is  possible,  ho^.jever ,  that 
mixing  by  drift  currents  accounts  for  the  larger  HLD  values  observed 
at  shallov7cr  depths,  i.e.  for  wind  speeds  less  than  12  kuccS/  r  so. 

The  modt:ls  of  wave  mixing  should  be  considered  as  indices  of  mixing 
efficiency,  not  as  precise  physical  descriptions  of  the  vertical  mixing 
process,  do  one  has  been  able  to  describe  the  exact  process,  but  the 
vjave  mixing  approach  a.ppears  to  be  the  most  promising  thus  far  advanced. 

It  is  difficult  to  say  from  one  test  which  model  pf  \;ave  mixing  as 
presently  v/ritten,  if  any,  V7ill  prove  to  be  useful  in  a  forecasting  scheme 
Laevastu*s  equations  give  the  closest  appro::imation  to  an  upper  limit  on 
the  observed  mixing.  Neumann's  energy  ratio  fits  the  observed  distribu¬ 
tion  better  chan  the  other  equations,  hov/ever  Geary's  model  takes  into 


account  thcrnocline  stability,  a  necessary  paraiiicte^  ii  other  couoitions 
aad  locations  arc  ever  to  be  investigated. 

The  time  lag  between  observed  \7ind  and  resulting  niixing  is  a 
provident  phenomenon  from  a  forecasting  standpoint,  and  must  be 
considered  when  verifying  calculated  mixed-layer  depths. 

Hourly  and  other  short-term  randoia  fluctuations  of  MLD  about  a 
mean  value  prevent  the  forecasting  of  exact,  ins tantancous  mixed-layer 
depths,  and  vjill  have  to  be  considered  for  forecasting  and  investigative 
purposes . 

The  most  important  factor  requiring  investigation  is  the  problem 
of  convective  mixing*  At  present,  this  process  has  no  model  to  des¬ 
cribe  its  effects,  yet  it  appears  to  be  the  singular  mechanism  domina¬ 
ting  the  depth  to  which  mixing  takes  place  during  most  of  the  annual 
cycle.  The  use  of  a  high-speed  computer  for  a  running  computation  of 
the  heat  balance  in  the  surface  vzaters  may  prove  to  be  of  valuable 
assistance  in  describing  short-term  fluctuations  of  mixed-layer  depth. 
Future  research  should  also  be  applied  to  improving  and  testing  models 
of  wind  mixing. 
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Figure  1 


Tims  Series  of  V/ind  and  MLD, 
Ocean  Station  Pane,  Juna,  1958 


June,  1958 


I 

20 


Wind  (knots) 


C 


29 


Wind  and  t-ff.D,  Ocean  Station  Pana,  June  through  St^ntcmi'yer,  1958 
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Figure  2 
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The  Effects  of  Convection  on  MT>D 
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APPLNDIX  I 


gi::a.iy  s  model 

I'icchaiilcal  ui::ip»7  ol  water  particles  eacross  the  tacrmoclinc  must 
iaaibited  when  the  vertical  forces  of  particle  laotion  are  opposed 
bl  t'lc  buoyant  force  of  the  denser  water  belov"  the  tnermoclinc.  Based 
upon  tais  consideration  it  is  possible  to  establish  a  simple  model  for 
mccheuical  mixing  by  v;avc  action  vjith  the  follow^ing  assumptions: 

1)  A  homogeneous  la\‘ci  of  V7ater  wata  density  yO  (the 
mixed  la)'er)  overlies  homogeneous  water  with  a  greater 
density,  JD  .  The  difference  in  density  between  the 

Vjo  layer-  is  caused  'solely  by  their  temperature  difference. 

2)  Subsurface  particle  motion  is  in  accordance  w’ith  simple 
Airy  \7ave  theory  for  deep  V7ater,  i.e,  the  particles  move 
in  circular  orbits  with  their  diameters  decreasing 
exponentially  \7ith  increasing  depth. 

3)  Penetration  of  the  thermocline  interface  stops  when  the 
vertical  forces  of  a  particle  in  motion  are  exactly  balanced 
b\  t:he  buoy^.nt  force  The  balance  of  forces  in  such  a  system 
are  given  by: 


^3 


Aa-  s  =  p's 


(2) 


2  - 

uAiere yOQ  is  the  gravity  force,  yO ^  is  tlie 

maximum  force  due  to  the  vertical  component  of  acceleration 

/ 

of  the  particle  in  orbital  motion,  and  yO  ^  is  the  buoyant 
force . 

The  force  balance  of  (2)  is  correct,  but  it  must  be  understood 
that  the  acceleration  term  0~  ^  applies  to  laminar 
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conditions  of  wave  noti  and  is  used  here  only  as  an  approximation 
to  the  magnitude  of  vertical  particle  acceleration  in  what  must  be 
actually  turbulent  moticu.  The  laminar  acceleration  term  is  utilized 
in  an  attempt  to  resolve  in  a  general  way  the  difficulties  of  this 
simple,  mechanistic  appr'^ach. 

Equation  (2'^  can  be  evaluated  implicitly  for  depth  as  a  function 
of  wind  speed  alone  if  the  stability  of  the  thermocline  is  assumed 
constant,  i.e.,  ^ P  =  constant^ for  all  depths  and  the  various 

surface  wave  character  i-:tic^  are  calculated  according  to  Neumann's 
spectrum  of  wave  generation  by  wind. 

t 

The  stability  term,  ,  was  evaluated  from  the  data  at  ocean 

P 

station  Papa  for  June  and  July  of  1958,  and  \7as  found  to  be  approx¬ 
imately  between  the  liiric^;  of  0.0001  and  0.0004.  Using  EO  to 
obtain  the  wave  character ii tics  T”  ,  ,  and  A/  for  fully  developed 

seas,  table  1  gives  the  computed  values  of  HLD  for  various  vjindspeeds 
and  two  stability  valuer?..  Expressing  (  2  )  as 

^  ~  act""  ^  P  J  ■>  (3) 

a  sample  calculation  follows:  for  a  given  windspeed  7”,  Z- , /V are  speci¬ 
fied  as  above.  A  s  sume  PZP  -  0.0004,  A  -  Jd.  j  CT=  RIT  ,  j 

^  ’2  ^  2: 

solve  (  3  )  for  ,  which  represents  the  MLD. 
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Table  1 


Geary  s  MLD  Values 

l.’iricls peed  (knots)  ...ed-Layer  Depta  (iacters) 

=  0.0001  ■  O.GOO^f 




10 

8.C 

6.9 

12 

12.8 

10.0 

\h 

17.6 

14.0 

16 

23.3 

18.4 

10 

29.7 

23.6 

20 

36.8 

29.4 

22 

45.0 

36.0 

Ih 

54.0 

43.0 

26 

63.5 

51.0 

23 

73.7 

59.4 
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APPENDIX  II 


LAEVASTU’S  MODEL 


Laevastu  give?  a  tentative  formula  for  the  depth  of  mixing 

by  waves;  uufor tuna^ el v  he  does  not  give  the  derivation,  but  states 

that  he  arrives  at  ^he  expression  by  using, 

the  relations  for  trochoid  waves  for  computing  the  velocities 
of  water  particles  at  various  wave  height  and  at  various 
depths,  the  depths  ^i^ren  by  Neumann  J  ?  where  the  total 
V7ave  energy  ha?  decreased  to  five  percent  of  its  value  at  the 
sea  surface  and  a^sa.ming  that  the  mixing  by  waves  is  negligible 
at  approximately  a  depth  where  the  diameter  of  the  orbital 
paths  is  smaller  than  10  centimeters. 

The  formula  is 


/^.S  /Vjr  . 


Laevastu  uses  the  foIlcv7ing  equation  for  computing  significant 
wave  height: 


_  0-0003  [^0+ C'^~  ~^)J 

^  /=^  /  (5) 

It  is  interesting  to  note  that  (5)  incorporates  the  sea-air  tempera¬ 
ture  difference  parameter,  while  other  formulas  for  significant  wave 
height  do  not.  Depth  of  the  mixed  layer  calculated  from  (4)  and  (5) 
is  given  in  table  2.  assuming  =  0,  F  =  100  kilometers,  and 

t  -  24  hours, 
a 

MLD  obtained  ti om  formula  (4)  is  greatly  influenced  by  the  equa¬ 
tion  used  to  calculate  significant  wave  height.  As  an  illustration  of 
this  point,  table  .3  give-  MLD  fcr  various  vjind  speeds  using  formula  (4), 
but  \7iLh  significant  wave  height  computed  according  to  Neumann’s  equation. 


26 


Table  2 


La^vastu  s  ITLD  Vj  lues 


\viud  speed 

hixcci-Layci:  Depth 

rii«3 ti 3  see 

kao  l:  s 

i.ic  tors 

3 

5.8 

3.9 

/j 

7.8 

6.6 

5 

9.7 

9.9 

7 

’3.6 

17.6 

8 

15.5 

22.3 

9 

17.5 

26.3 

10 

19.^. 

31 .8 

12 

23.3 

42 . 3 

Table  3 

Modified  LaevasCu  MLD  Values 

Viad  s/eed 

Mixed-Layer  Deytli 

luCtC^'S  SGC"*^ 

knots 

uetcr s 

3 

5.3 

1.4 

4 

7.8 

2.8 

5 

9.7 

5.0 

■/ 

13.6 

1 1 . 4 

U 

15.5 

16.0 

J 

17.5 

21.5 

10 

19.  A 

17.9 

12 

23.3 

44 . 1 

2  7 


/*irL.;Di}.:  iii 
.T  ’i:yV.:.)'S  hOELL 


C.U  er.iejici,  . 


t.car.uian 


C^O  dfc’ IV-:.  .. 


enc-tj,  L  ,^pLh  i  f'  d 

\  cl.  S  ctS  t*  IJiICL-  I' 


^1;  :;pectrum  of  v;ave  ;jcn<-ration  by  \7inc’ , 
:quat'ou  for  the  ratio  of  average  potential 
u  potential  energy  of  the  surface 

ind  ^'peed  and  depth: 


-h  25  1  e  ~ 

J  j 


(6) 


it  is  possible  to  i*. the  vertical  displacement  of  a  particle  at 
cc.iy  an’  for  a  /ivei  •.;ind?pecd  by  substituting  the  follo\;ing 

i.  v-ir.Cioasaip&  il’Lo  *6. 


^  ^  5  (7) 

(8) 


c  uuann  con s  i  e er  £  t b u 


\7\je  aoLion  is  ner'i  i-^i 


u. 


is  eqv.al  to  .  0 


^,eunUT  . 


v.on  the  ratio  of  becomes  equal 
Tjo 

'^le.  Solving  (6)  for  the  condition 
5  table  4  is  obtained. 

Table  4 

i  Dep^h  of  Negligible  Nave  Motion 


to  0.05, 
:here 


'./iiid  soced 


meters  sec 


Mixed -Layer  Depth 
meter  s 


2 

0.42 

5 

5.7 

10 

22.9 

15 

52 

20 

92 

25 

143 
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